Open quantum random walks: bi-stability on pure states 
and ballistically induced diffusion. 
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Open quantum random walks (OQRWs) deal with quantum random motions on the line for 
systems with internal and orbital degrees of freedom. The internal system behaves as a quantum 
random gyroscope coding for the direction of the orbital moves. We reveal the existence of a 
transition, depending on OQRW moduli, in the internal system behaviours from simple oscillations 
to random flips between two unstable pure states. This induces a transition in the orbital motions 
from usual diffusion to ballistically induced diffusion with large mean free path and large effective 
diffusion constant at large time. We also show that mixed states of the internal system are converted 
into random pure states during the process. We touch upon possible experimental realisations. 
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Introduction. Random walks j3] are ubiquitous in our 
understanding of physical phenomena with plethora of 
applications in biology or economics. They are instru- 
mental in mathematics and computer science. Quantum 
generalisations have been considered decades ago [I] and 
they find numerous applications in quantum computa- 
tion or quantum cryptography [5]. They have recently 
been experimentally implemented [6]. Drastically influ- 
enced by quantum interferences, quantum random walks 
behave very differently from their classical analogues, for 
instance they do not diffuse in the same way. 

OQRWs were introduced [7J using concepts from dy- 
namical maps [H] aiming at incorporating decoherence ef- 
fects. They specify random motions of quantum systems 
with both internal and orbital degrees of freedom (d.o.f), 
and these moves depend on interactions with quantum 
coins. Contrary to quantum random walks, OQRWs im- 
plement resettings of the quantum coins at each time 
step, and this difference has profound consequences. 

Studying classes of OQRWs we find a transition in 
their behaviours separating usual diffusions from ballis- 
tically induced diffusions with large mean free path be- 
tween trajectory flips. Of course diffusion is always due to 
ballistic behaviours at a small enough scale, what matters 
is the time separation between flips. In OQRWs, these 
are not due to disordered collisions but to abrupt tilts 
of the internal gyroscope induced by the interaction with 
quantum coins and their measurements. Behaviours in 
the ballistic regime are consequences of random switches 
of the internal state between unstable pure states. 

Open quantum random walks. To be closer to possible 
experimental realisations and to quantum trajectory the- 
ories [21 [TO], we define OQRWs using a picture slightly 
different but equivalent to [7] in which the system in- 
teracts recursively with identical quantum coins, called 
probes [11]. The probe Hilbert space rl p is choosen to 
be two dimensional with a specified basis {|±) p }. At 



each time step, the system interacts quantum mechan- 
ically with one sample of identically prepared copies of 
the probe on which a measurement is performed after the 
interaction period. The system-probe interaction is such 
that if the out-going probe is measured in the state \+) p 
(resp. \—) p ) the system moves by one step to the right 
(resp. to the left) along the line, and this move is ac- 
companied by a modification of the internal d.o.f 's. The 
system position is slave to the measurement out-puts. 

Although experimental realisations of OQRWs do not 
yet exist we may contemplate possible scenarios in which 
the orbital states |n) refer either to energy levels, num- 
bers of photons or localised positions. One may imagine 
using ions trapped in harmonic potentials, as in [121 113) . 
each ion being possibly in two states with different angu- 
lar momenta, and photons as probes. For an appropri- 
ately adjusted frequency and linearly polarised in-going 
photons, the ion-photon interaction may induce internal 
flips and energy shifts conditioned on the measurements 
of out-going photons [14] . One may also imagine using 
cold atoms with internal d.o.f 's and localised on potential 
lattices, as in [5J, and probing them coherently with pho- 
tons |14) . If one is only interested in the internal system 
|15) . a set-up dealing with recursive couplings of a Qbit 
to series of probe Qbits, as in [TO], may be considered. 

To make this description concrete, let H S ys '■— Hc®Ho 
be the system Hilbert space, with rl c and rl respectively 
associated to the internal and orbital d.o.f 's. We take rl c 
finite dimensional and T~L ~ C z with orthonormal basis 
{|n)o, n 6 Z}. Let U be the unitary operator acting on 
H S y S <8> 'Hp coding for the system-probe interaction. We 
demand that its action on states \ip) c <& \n) <8 \4>)p gives 
the entangled normalised states 

(B + \i/>) c ) <E)\n+ 1) ® \+) p + (B_ |V>) C ) ® \n - 1) ® h) p , 

for any \tp) c eH c . Unitarity imposes b\_b + +bIb- = i. 
OQRWs consist in iterating system-probe interactions 



2 



and out-going probe measurements. Since the latter are 
random with probabilities governed by quantum mechan- 
ics, this yields stochastic evolutions called quantum tra- 
jectories [HI [TO]- If the system density matrix is initially 
localised in the orbital space, say po <X> |£o)o(^o|, it re- 
mains so after each iteration with internal density matrix 
p n and orbital position x n . These are randomly updated, 



Pn ® \x n ) o(%r, 



B±p n B\ 



\x n ± l)o(x n ± 1 



(1) 



with probability p^ ■= tr-n c (B±p n B±). The process 
n — > (p n , x n ) is Markovian on a probability space whose 
events are the recursive out-put probe measurements. By 
construction the mean system density matrix evolves ac- 
cording to the OQRW quantum dynamical map [7j, and 
the mean internal density matrix p n := E[/0„] satisfies 
Pn+i = B + p n B\ + B-p n B^_. In absence of internal 
d.o.f 's OQRW behaviours parallel those of classical ran- 
dom walks. We take T-L c ~ C 2 and represent the internal 
system by an effective spin one-half gyroscope. 

Heuristics. Solutions of the unitary constraint may be 
parametrised as B± — U±M± with U± unitary and M± 
hermitian with M 2 + Af 2 = I. These are the OQRW 
moduli. In the numerical simulations we look for OQRWs 
generated by B+ = and = J -1 ( ) 

with S = V u 2 + v 2 + r 2 + s 2 , and in the scaling limit by 
B± = ^[1 ± ^/eN + e(iH ± ± M - §iVtjV) + o(e)] with 
e a small parameter and H±, M hermitian but not N. 
We take H := \{H+ + H-) = u Q a 2 and N = a a 3 with 
0-1,2,3 -j- ne usua \ Pauli matrices. 

Numerical simulations reveal the existence of differ- 
ent regimes for OQRWs corresponding to a 2 /oj a below 
or above a critical value. For a 2 /loq small enough, the 
position x n is nearly brownian and the internal density 
matrix p n oscillates almost regularly, see Figj2j More 
interesting behaviours occur for a 2 /luq above the critical 
value, see Fig[T] The position x n follows a random seesaw 
trajectory, with tiny fluctuations, whose slopes are deter- 
mined by the internal state which fluctuates around two 
unstable fixed points and toggles randomly from one to 
the other. The abrupt changes in the position moves are 
due to the random flips of the internal gyroscope. The pa- 
rameter a 2 /luq controls the mean free path between flips. 
Although ballistic on this time scale, the position is diffu- 
sive on larger time scale. Whatever the initial value, the 
internal density matrix converges rapidly to pure states, 
so that the fixed points are also pure states. It is quite 
remarkable that series of indirect probe measurements 
project mixed states on pure states. The progressive col- 
lapses elegantly observed in [T7|, and proved in [TH], is a 
particular illustration of this phenomena, but in OQRW 
context the target states keep on evolving randomly. 

This peculiar behaviour bears similarities with that 
of a noisy particle in a double well potential subject to 
Kramer's transitions from one well to the other. This is 
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FIG. 1: OQRW generated by B± as in the text with u = 1.1, 
v = 1.00 and r = -s = 0.00015: (a) position X t , (b)k,(c) 
a 1 & a 3 -components of p t , (d) determinant At. 



the picture that we are going to make explicit in the fol- 
lowing. The convergence towards pure states can be un- 
derstood as follows. Let A„ := det p n be the determinant 
of the internal density matrix. It vanishes for pure states. 

,1/2 



1 /2 

A simple computation shows that E[A„ ] = 
c := det ^ (Ml) + det 5 (Mi) < 1, and Iim n _, 



c" A ' with 

ooE[Ay 2 ] = 



1/2 

0. Actually we can prove that linin^oo A n = almost 
surely using the sub-martingale convergence theorem of 
probability theory |19j . Indeed, computing the mean of 

1/2 

A n+1 conditioned on the n-first out-put measurements 



gives E[A£ft|.Fn] = cA n ,z < A 
sub-martingale, and since it is bounded, it converges al- 
most surely and in L 1 . The limit is zero and the internal 
density matrix localises on pure states. As numerically 
observed, the convergence is very fast. 

In a continuous limit, the mean system density matrix 
reads Jdx p(x,t) ® |ar) with p(x,t) := tr-n c p(x,t) the 
probability density to find the system at position x at 
time t, and p t '■= fdx p(x,t) the mean internal state. At 
each time step dt, it is updated using OQRW rules Q, 

p(x, t + dt) = B-p(x + dx, t)B[_ + B + p(x - dx, t)B\_. 

A continuous limit exists if one imposes the scaling rela- 
tion e = dt = dx 2 |20j . Taylor expansion then gives: 

dtp - \dl P - (Nd xP + d xP N^) + i[H, p] + L N (p), (2) 

with Lindbladian L N {p) := pN - \{NN^p + pNN^). 

mixes pieces from diffusive Fokker-Planck equa- 
tion and from Lindbladian quantum evolution for p t |21j . 
The term (Nd x p + d x pN^) is at the origin of the ballistic 
behaviour seen in FigJT] and of the large effective diffu- 
sion constant but the hamiltonian term is required for the 
tilting effect. The probability p(x, t) is not associated to 
a Markov process and does not satisfy a linear equation 
but it becomes gaussian at large t. 
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so that aI/ 2 is a 
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FIG. 3: A yt -trajectory for a 2 > Uq. 



FIG. 2: OQRW generated by B± as in the text with u = 1.005, 
v = 1.00 and r = -s = 0.00015: (a) position X t , (b)fc(c) 
a 1 & cr 3 -components of p t , (d) determinant A t . 



Scaling limit. Let us now make precise the heuristic 
description by deriving the stochastic differential equa- 
tions (SDEs) governing OQRWs in the scaling limit. 
OQRWs are defined on the probability space whose 
events are the series (si,S2,---) with Sk = ± depend- 
ing whether the k-th out-going probe is measured in the 
state |±)p. Functions which depend only on the n first 
, s n ) define a natural filtration T n [19] . and 
+1= ±}\J r n] = tr(S±p„_Bj_) are the probabili- 



data (si, • • 

pt ■■= nhs 

ties for s n +i = ± conditioned on the value of the internal 
state at the n-th step. A quick and neat way to obtain 
the scaling limit consists in decomposing the process p n 
as a sum of a martingale M n plus a predictable process 
O n - This is called a Doob decomposition [19]. In the 
scaling limit the martingale (resp. predictable) contribu- 
tion converges to the noisy source (resp. the drift) of the 
SDEs. Eq.0 may be tautologically written as 



Pn + l = Pn +) l{ Sn+1 =+} + P { n 5 I{ Sn+1 =-}, 



,(±) 



and x n +i - x n = I{ s „ +1 =+} - I{ Sn+1 =_}, with p„ 
B±p n B±/pn ■ By construction the Doob martingale is 
M n = YX=\ Wk witn nk := P k ~ ^[Pk\Fk-i] given by: 



2^ = {Pk ~ Pi (hs k+1=+} -Pt+Pk- hs k+1 =- } ) ■ 

The predictable process is defined by complementarity 
O n '■= Pn — Mn- Taking the scaling limit e — > 0, 
t = ne fixed, is a matter of Taylor expanding dM t := 
M n+1 -M„, dp t := Pn+i-Pn and dX t := y/e(x n+ i-Xn)- 
Identifying e with dt, we get dM t = Z)jv(pt) dB t , and 



dp t = (i[H,pt] + L N (p t ))dt 
dX t = U N { Pt )dt+ dB t , 



D N {p t )dB u (3) 
(4) 



with B t a normalised Brownian motion. 
pATf - p u N (p) and U N (p t ) := tr(JVp 4 
prisingly, eq.(|3| is of Belavkin's type [2 



Z?jv(p) :=iVp + 
piV t ). Not sur- 
31531. The drift 



in eq.Q is governed by the internal state and this is re- 
sponsible for the behaviours observed in Figs(T^[2] 

Bi-stability and ballistic diffusion. We take H = ujq a 2 
and N — a a . Eqs.(|3|4[) are then compatible with reality 
of the internal density matrix. We para metri se it as pt = 



^(I + qia 1 +q 3 a 3 ) with ql + q 2 < 1. Eqs.l3 4|) then reads: 



dq 3 = 2wq q\ dt + 2a(l — q 2 ) dB t 

dqi = -2((J q 3 + a 2 q x ) dt - 2a qiq 3 dB t 

To check convergence to pure states, let A t := det p t . We 
have dA t — — 4A t [a 2 (l— q 2 )dt+aq 3 dB t ] with non positive 
drift so that A t is a sub-martingale |19j . It converges 
quickly to 0, so we may describe p t as a pure state, q± = 
sin 8, q 3 = cos 9. The angle 9 satisfies 



d9 t 



-2(coo + o 2 sin 9 t cos 9 t )dt — 2a sin 9t dB t (5) 



The behaviour of 9t is quantitatively different depend- 
ing whether a 2 ^ cj , and this corresponds to the two 
regimes we mentioned. For a 2 < luq, 9 t rotates randomly 
but regularly enough around the unit circle, so that the 
internal state pt oscillates almost regularly. For a 2 > uiq, 
9 t is trapped during random periods in the vicinity of 
6*_ ~ 0~ or 9* + ~ 7T _ . The points 6± are the minima 
of the effective potential obtained from eq.([5| once cor- 
rectly normalised. Although it fluctuates, 9 t turns pre- 
dominantly clockwise (for uj > 0) around the unit circle, 
never crossing back or n anticlockwise. 

To make this description quantitative, let y t := 
— log | tan#t/2|. It satisfies a normalised SDE with con- 
stant noise source, dyt = 2a dB t — V \y t )dt with potential 

V{y) = — 2(±wo sinhy + 2a 2 log cosh y). 

The above sign is that of tan# t /2, i.e. +/— for 9 t on 
the upper/lower half unit circle. What happens in these 
two sectors is symmetrical, so we concentrate on the up- 
per sector. The potential shape is that of a cubic like 
function but it is exponentially large for large |y|, i.e. 
V(y) ~ — u>o sign(y)e' J '' . It posseses a minimum and a 
maximum for a 2 > ojq, and none if a 2 < u>q. The mini- 
mum is at y* + ~ — 2a 2 /wo for large a, i.e. tan#^j_ ~ e~ v + 
so that 9+ is close to ir~, with I^nin — — 4a 2 log a 2 /wq 
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and Knax — 0. When 9 t enters the upper sector, it does 
it from 7r. For y t this corresponds to — oo, so that yt ex- 
periences an exponentially steep down ramp that it can 
never climb back, and this means that 6 t never escapes 
the upper sector from tt but only from 0. Going down on 
the ramp, y t reaches the potential minimum and spends 
time fluctuating around there, and this means that 8t 
fluctuates around 9* + . At a random time Tfli p , large fluc- 
tuations allow y t to cross the energy barrier in a Kramer's 
like process. Once this has happened, y t is again on a 
steep ramp that it steps down to +oo, and this trans- 
lates to 9t moving toward + and crossing it irreversibly 
towards the lower sector. The process then starts on the 
lower half circle and repeats itself. See Fig(3l We esti- 



J.M. Raimond for discussions, especially on possible ex- 
perimental scenarios. 



mate the mean flip time as E[Tflj p ] 



,AV/i 



by 



Kramer's rule, and a more precise study allows us deter- 
mine the probability distribution of Tfli p /a 2 . 

The internal state drives the system position via eq. Q 
which reads dX t = 2a cos 0t dt + dB t . The slopes of the 
seesaw profiles of X t are 2a cos 6^ — =F2a. Fluctuations 
are negligible for large a but the noise is instrumental 
for tilting from one slope to the other via Kramer's tran- 
sitions. The mean system density matrix heuristically 
introduced above is rigourously defined by 



dxp{x,t) ® \x) (x\ := E[p t ® \X t ) (X t 



(6) 



Routine applications of stochastic ltd calculus [TH] show 
that the SDEs ([30) imply eq.@ for p(x,t). Eq.Q is of 
Lindblad form on H c <8 L 2 (M). It may be used to check 
that Xt/y/t becomes gaussian at large time, in a way 
compatible with the central limit theorem of J2JJ 023] , and 
E[AT t 2 ] = D c g t with effective diffusion constant D c g = 
1 + 4a 4 /cjq. The factor I comes from the bare diffusion 
constant |20j while the second term, which dominates for 
large a, is induced by the ballistic seesaws. 

Conclusion. The transition from usual diffusion to bal- 
listically induced diffusion is an echo of the internal gy- 
roscope behaviours. In the ballistic regime the internal 
state switches randomly between two pure states in a 
way similar to Kramer's transition. Since the system po- 
sition is slave to the out-put measurements, our results 
about convergence from a mixed state to pure states and 
about random flips between them apply to the coupled 
probe plus internal spin system without considering or- 
bital d.o.f's. More details will be given elsewhere [2"rJ] . 
In the ballistic regime the effective diffusion constant is 
much larger than the bare one, and one may wonder 
about other scenarios of ballistically induced diffusion 
providing large effective diffusion constants. Finally, one 
may muse about possible applications of the rapid and 
irreversible internal state flips to switch devices [27]. As 
such these switches are random but one may imagine 
changing on demand the height of the energy barrier. 
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